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Abstract 

Motivated by the work of Anstee, Griggs, and Sali on forbidden submatrices 
and the extremal sat-function for graphs, we introduce sat- type problems for 
matrices. Let J-" he a family of /c-row matrices. A matrix M is called J-'- 
admissible if M contains no submatrix F G (as a row and column permutation 
of F). Moreover, M is -saturated if M is ^-admissible but the addition of any 
column not present in M violates this property. In this paper we consider the 
function sat(n, J^) which is the minimum number of columns of an .7^-saturated 
matrix with n rows. We establish the estimate sa.t{n,J-) = 0{n^~^) for any 
family of /c-row matrices and also compute the sat-function for a few small 
forbidden matrices. 



1 Introduction 

First, we introduce some notation. Let the shortcut 'an m x n-matrix' M mean a 
matrix with m rows (which we view as horizontal arrays) and n 'vertical' columns. 
Usually, we restrict entries to only two values, and 1; in all other cases, we indicate 
the range explicitly. For an n x m-matrix M, its order v{M) = n is the number of 
rows and its size e{M) = mis the number of columns. We distinguish the expressions 
like 'an n-row matrix' and 'an n-row' standing respectively for a matrix with n rows 
and for a row containing n elements. 

For an n X m-matrix M and sets A C [n] and B C [m], M{A, B) is the |^| x \B\- 
submatrix of M formed by the rows indexed by A and the columns indexed by B. 
We use the following obvious shorthands: M{A, ) = M{A, [m]), M{A, i) = M{A, {i}), 
etc. For example, the rows and the columns of M are denoted by M(l, ),..., M{n, ) 
and M{,1), . . . , M{,m) respectively while individual entries - by M{i,j), i G [n], 
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A matrix F is a suhmatrix of a matrix A (denoted as F C A) if deleting some 
set of rows and columns of A we can obtain a submatrix which is a row/column 
permutation of F. The transpose of M is denoted by M"^ (which we use mostly to 
denote vertical columns for typographical reasons); (a)* is the (horizontal) sequence 
containing element a i times. The n x (mi + m2)-matrix [Mi,M2] is obtained by 
concatenating an n x mi-matrix Mi and an tt, x Tn2-matrix M2. We write N = M 
to say that is a column/row permutation of M. Thus, C M if = M{A,B) 
for some index sets A and B. The complement 1 — M of a matrix M is obtained by 
interchanging ones and zeros in M. The characteristic function xy of y C [n] is the 
n-column with ith entry being 1 ii i and otherwise. 

Many interesting and important properties of matrices (that is, classes of matrices) 
can be defined by listing forbidden submatrices. (Some authors use the term 'forbid- 
den configurations' then.) More precisely, given a family JF of matrices (referred to 
as forbidden) , we say that a matrix M is J-" -admissible (or T-free) if M contains no 
F G as a submatrix. A simple matrix M (that is, a matrix without repeated 
columns) is called J^-saturated or (J^-critical) if M is jF-free but the addition of any 
column not present in M violates this property; this is denoted by M G SAT(r2,jF), 
n = v{M). Note that, although the definition requires that M is simple, we allow 
multiple columns in matrices belonging to JF. 

One well-known extremal problem is to consider forb(n,jF), the maximum size 
of a simple jF-free matrix with n rows or, equivalently, the maximal size of M G 
SAT(n, JF). Many different results on the topic have been obtained, some of them will 
be mentioned in the due course, but we do not even try to give a comprehensive survey 
here. We would only mention a remarkable fact that one of the first forb-type results, 
namely formula ([T]) here, proved independently by Vapnik and Chervonenkis |21j . 
Perles and Shelah [12], and Sauer [IB], was motivated by such different topics as 
probability, logic and a problem of Erdos on infinite set systems. 

The forb-problem is reminiscent of the Turan function ex(n,jF): given a family 
JF of forbidden graphs, find the maximum size of an ^-free graph on n vertices 
not containing any member of JF as a subgraph (see e.g. surveys [13 [2D])- Erdos, 
Hajnal, and Moon [TT] considered the 'dual' function sat(?7,,jF), the minimum size 
of a maximal ^-free graph on n vertices. This is an active area of extremal graph 
theory, see [6l Section 3] for a short overview. 

Here we consider the 'dual' of the forb-problem for matrices. Namely, we are 
interested in the value of sat(n, JF), the minimum size of an jF-saturated matrix with 
n rows: 

sat(n, J^) = min{e(M) : M G SAT(n, J")}. 

We decided to use the same notation as for its graph counterpart. This should not 
cause any confusion as this paper will deal with matrices alone. 

Obviously, sat(n,^) < forb(n,jF). If JF = {F} consists of a single forbidden 
matrix F then we write SAT(n, F) = SAT(n, {F}), etc. 

A simple n-row matrix M is monotonically J-" -saturated, denoted M G m-SAT(n, JF), 
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if the addition of any new n-column C to M creates a new forbidden submatrix, 
that is, there is yl C [n] such that M{A^ ) is jF-free while [M, C](v4, ) is not. Clearly, 
m-SAT(n, JF) D SAT(n, T), so m-sat(n, JF) < sat(n, where as always m-sat(?7,, JF) = 
min{e(M) : M G m-SAT(n, JF)}. This is the natural analog of the m-sat-function for 
graphs, which was studied already in [TT]. It is useful in proving lower bounds on 
sat(n,jF) via induction on n. 

There is an obvious generalization of these problems when we consider the class of 
[0, l]-matrices (matrices whose entries can assume (/ + l)-values from {0, ...,/}) with 
the above definitions going practically unchanged. In this case we will use symbols 
like SAT(ri, JF; etc., whereas the default SAT(n,jF) = SAT(n,jF;l) is the usual 
notion. 

By we denote the simple k x -matrix consisting of all A;-columns with exactly 
/ ones and by i^'^-the k x 2^ matrix of all possible columns of size k. Naturally, T^r' 
denotes the k x f[k, /)-matrix consisting of all distinct columns with at most / ones, 
etc, where we use the shortcut 

/(..)^(^). (")... ..(:). 

Vapnik and Chervonenkis [H], Perles and Shelah [T^, and Sauer [TS] showed 
independently that 

ioTh{n,K,) = f{n,k-l). (1) 

As formula ([1]) plays an important role in our consideration of the sat-function, 
we give a proof, for the sake of completeness, of a more general result which is of 
independent interest. Namely, in the class of [0, /]-matrices, we compute the maximum 
size of a X[-free n-row matrix, where K^. is the /c-row matrix of size + consisting 
of all distinct A;-columns made of symbols from [0, 1] = {0, ...,/}. Our proof is based 
on the ideas of Frankl introduced in [13] . 

Theorem 1.1 For any I > 1, k > 1 and n > k — 1, 

iorh{n,Ki;l)=J2l''~'('')- (2) 

Proof. To prove the lower bound, consider the matrix made of all columns containing I 
at most k — 1 times which is obviously X^-admissible. 

It remains to prove the upper bound. Consider any i^'^-admissible matrix M of 
order n and size m. Fix an integer i, 1 < i < n. Divide the set of columns of M 
into disjoint sets Ci, . . . ,Cp, 1 < p < m, such that all columns in Cj, I < j < p, agree 
except in the ith row. Clearly, \Cj\ < I + 1. Let be obtained from Cj by replacing 
the set of entries corresponding to the ith row in columns of Cj by {0, . . . , \Cj\ — 1}. 
Let Ml be a matrix consisting of all columns from Cj over all 1 < j < p. 

Claim. If Ml D Kl then M D Kl 
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Suppose Ml D . Let R be the subset of k rows of such that M/(i?, ) 3 
If z is not in R then M also contains (since M and are identical apart from 
row i). Thus, we may assume that i is in R. Let C be any column of . Let be 
the same column with I in the zth entry. Then some column of M- gives when 
restricted to R. And has an / in the zth entry, so if is in is also in M, 

and every column obtained from by putting anything in the entry is also in 
M. But then C is a column of M when restricted to R. Thus M D X^, as required. 

In order to prove the theorem we recursively repeat the claim for every row i, 

1 < < n, obtaining a matrix M' of size m. Observe that if any column of M' has k 
entries equal to / then M' ^ K\. and so by the claim M D a contradiction. Hence, 
every column of M' has at most k — 1 entries equal to / which gives the theorem. □ 

2 General Results 

Here we present some results dealing with sat(n,jF) for a general family T. 

The following simple observation may be useful in tackling these problems. Let 
M' be obtained from M G SAT(ri, T) by duplicating the nth row of M , that is, we let 
M'([n], ) = M and M'(n+ 1, ) = M (n, ). Suppose that M' is J^-admissible. Complete 
M', in an arbitrary way, to an jF-saturated matrix. Let C be any added {n + 1)- 
column. AsbothM'([n],) and M'([n-l]U{?i + l}, ) are equal to M G SAT(n, J"), we 
conclude that both C([n]) and C([n— 1] U {?7,+ 1}) must be columns of M. As C is not 
an M'-column, C = (C, b,l — b) for some {n — l)-column C and b G {0, 1} such that 
both (C, 0) and (C, 1) are columns of M. This implies that sat(ra+l, JF) < e{M)+2d, 
where d is the number of pairs of equal columns in M after we delete the nth row. 

The above argument works for the m-sat-function. Namely, duplicate a row of 
some M G m-SAT(n, JF) and add one by one missing columns so that no new forbidden 
submatrix appears; clearly the resulting matrix is monotonically jF-saturated and, like 
above, we add at most 2d extra rows. In particular, the following theorem follows. 

Theorem 2.1 Either m-sat(n, JF) is constant for large n or m-sat(n, JF) > n + 1 for 
every n. The analogous statement is true for the sat-function if no matrix in JF has 
two equal rows. 

Proof. If we have some M G m-SAT(n, JF) with at most n columns then a well-known 
theorem of Bondy [7] (see, e.g., Theorem 2.1 in [5]) implies that there is i G [n] such 
that the removal of the ith row does not create two equal columns. Hence, the 
duplication of the ith row gives a monotonically .F-saturated matrix, which implies 
m-sat(n+l, JF) < m-sat(n, JF) and the first part follows. The same argument establish 
the second part as the duplication of a row cannot create a forbidden submatrix by 
the condition on JF. □ 

Suppose that JF consists of k-Tow matrices. Is there any good general upper bound 
on forb(n,jF) or sat(n,jF)? There were different papers dealing with general upper 
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bounds on forb(n, JF), e.g., by Anstee and Fiiredi [1], by Frankl, Fiiredi and Pach [13] 
and by Anstee [1], until the conjecture of Anstee and Fiiredi [1] that forb(m, JF) = 
0{n^) for any fixed T was elegantly proved by Fiiredi (see [2] for a proof). 

On the other hand, we can show that sat(?T,, T) = 0{n''~^) for any family JF of k- 
row matrices (including infinite families). This is the matrix analog of the main result 
in [T7]. Note we cannot decrease the exponent of A; — 1 with the estimate remaining 
true for any JF; for example, sat(?7,, T^) = f{n, k — 1). 

Theorem 2.2 For any family of k -row matrices, sat(n, JF) = 0(n^~^). 

Proof. We may assume that is jF-admissible {i. e. every matrix of JF contains a pair 
of equal columns) for otherwise we are home by ([I]) as then sat(n, JF) < forb(n, Kk) = 
0(n'=-i). 

Let / G [0, A;] be the smallest number such that there exists m for which [mT^-', T^'] 
is not jF-admissible. (Clearly, such / exists as T^^ = is the complete matrix.) Let 
d be the maximal integer such that [mT^\ dTl, T^^] is jF-admissible for any m. Note 
that 6? > 1 as T^, T^^] = [mTj^'', T^-'] cannot contain a forbidden submatrix by 

the choice of /. Choose minimal m such that [mT^'', {d-\-l)Tl, T^'] is not jF-admissible. 

Suppose first that / < k. Given n, let be the n-row matrix corresponding to 
the following set system: 

H= [j G :E,6y?/ = J (mod n)}. 

ie[d-i] 

Here X^''^ = {Y C X : \Y\ = i} denotes the set of all subsets of X of size i. 

Note that any A E [n]^'^ is covered by at most d—1 edges of H as there are at most 
d — 1 possibilities to choose i E [n]\A so that A U {i} G H: namely, i = j ~ XlaeA ^ 
(mod n) for j G [rf — 1]. 

On the other hand, the set H' of all /-subsets of [n] covered by fewer than d—1 
edges of H has size at most 2(d—l) "-j^) . Indeed, if A E H' then, for some j E [d— 1] 
and X E A, 2x = j — J2a£A\{x} ^ (mod n) so, once A \ {x} and j have been chosen, 
there are at most 2 choices for x. 

Call X E [n]^'^) bad if, for some A E X^^\ 

\{Y E H :¥ D A, Y r]{X\A) = flS}\ < d-2. (3) 

To obtain a bad A;-set X, we either complete some A G if' to any k-set or take any 
/-set A and let X D A intersect some ii-edge covering A. Therefore, the number of 
bad sets is at most 

^(-<,-\)(/-^(:)<-<.-:-^«("'-')- 

Let M' = [A,T^]. Clearly, 

M'(A,) C [{-)T<\dTiTi^']. for any X E [n]W. 
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Hence, M' cannot contain a forbidden submatrix. Now complete it to arbitrary 
M= [M',M"] e SAT(n,J^). 

Suppose that e(M") C>(n'=-i). Then, by ([1]), Kk = M"{X,Y) for some X,Y. 
Now, remove the columns corresponding to Y from M" and repeat the procedure as 
long as possible to obtain > 0{n^^^) column-disjoint copies of Kk in M". If some 
X E [n]'^'^^ appears more than d times then (because T^{X, ) D ^T^^ for n > m + k) 
M{X, ) D [mT<\ {d+l)Kk\ is not J'-admissible. Otherwise, Kk C M"{X, ) for some 
good {i.e., not bad) X e [n]^''); but then A^(X, ) ^ {d - 1)T^ and M(X, ) contains a 
forbidden matrix. This contradiction proves the required bound for / < k. 

Suppose that / = /(JF) equals k; the above argument does not work in this case 
because M' ^ is too large. Let A^^ denote the {d + s) x d-matrix made of X[s+d]\{i}, 
i E [s + 1, s + d], where d = d{T) is as above; clearly, each A^^ is jF-free. Also let C* 
be obtained from a column C by adding d zeros to the end; to obtain M*, we apply 
this operation to every column of a matrix M; this increases the number of rows by 
d. 

Let JF* consist of all matrices M such that [M*, N^(m)] is not JF-free; clearly, this 
property is not affected by a row/column permutation of M. 
Note that T = [mT<'',T^] is not JF*-admissible as 

[T*, Nkm,) = [T, dT^] D [mT<^ {d + l)T',] 

so /(JF*) < k and, by the above argument, we can find L G SAT(n — d,J^*) with 
0(ri'^~^) columns. Now, M' = [L*,Nn-d\ is JF-free so complete it to an arbitrary 
M G SAT(?7,, JF). Let C be any added column; as [M', C]([n— rf], ) is jF*-free (otherwise 
M contains a forbidden submatrix), C{[n — d]) either is a column of L or equals 
((1)"-'^)^. Hence, e(M) < 2'^(e(Mi) + 1) and the theorem follows. □ 

3 Forbidding Complete Matrices 

Let us investigate the value of sat(n,A'fc). (Recall that Kk is the k x 2'''-matrix 
consisting of all distinct /c-columns.) We are able to settle the cases k = 2 and /c = 3. 
We will use the following trivial lemma a couple of times. 

Lemma 3.1 Each row of any M G m-SAT{n, Kk) , n > k, contains at least I ones 
and at least I zeros, I = 2^~^ — 1. 

Proof. Suppose that the first row M(l,) has mo zeros followed by mi ones with 
fTiQ > and / > mil. 

For i G [mo], let Ci equal the ith column of M with the first entry replaced by 
\. Then M' = [M, Cj] cannot contain Kk because the first row of M' contains too 
few I's while Cj([2,n]) is a column of M([2,n],). Therefore, Ci must be a column of 
M and mo = m-i. 

But then M has at most 2^^ — 2 columns, which is a contradiction. □ 
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The following result sounds as a rather natural (and not difficult) question if 
reformulated in the terms of set systems but we have not been able to find it in the 
literature. 

Theorem 3.2 For n> 1, m-sat(ri, K2) = sat(n, K2) = forb(n, K2) = n + 1. 

Proof. Suppose that the statement is not true, that is, there exists a monotoni- 
cally X2-saturated matrix with its size not exceeding its order. By Theorem 12. ![ 
m-sat(n, K2) is eventually constant so we can find an nxm-matrix M G m-SAT(?T,, K2) 
having two equal rows for some n G N. 

As we are free to complement and permute rows, we may assume that, for some 
i > 2, M(l,) = ... = M{i,) while M{j,) ^ M(l,) and M{j,) ^ 1 - M(l,) for 
any j G [i + 1, n]. Note that z < n as we do not allow multiple columns in M (and 
m > v{K2) — 1 = 3). By permuting the rows in M{[i + 1, n], ), we may assume that 
for some / G [z, n], we have M({1, j}, ) ^ if and only if j G [i + 1, /]. 

Let j G [i + l,?^]. By Lemma 13. the jth row M{j,) has both entries and 
1. By the definition of i, M{j,) is not equal to M(l,) nor to 1 — M(l,). It easily 
follows that there are fj,gj G [m] with M{ljj) = M{l,gj), M{j,fj) ^ M{j,gj); 
furthermore, we can find hj G [m] with M{l,hj) = 1 — M(l,/j). Let bj = M[j,hj)] 
we may assume that M{j, Qj) = hj. Furthermore, for j G [i + 1, /], we choose fj,gj, hj 
so that bj = M{j, 1), which is possible as M({1, j}, ) ^ K2. 

Now, as M G SAT(r7,, K2), the addition of the column 

c = ii,ioy-\k+,,...,bnf 

(which is not in M because C(l) 7^ C{2)) must create a new X2-submatrix, say in 
the xth and yth rows, some 1 < x < y < m. Clearly, {x, y} ^ [i] because each 
column of M([z], [m]) is either ((0)*)"'" or ((1)*)^. Also, it is impossible that x G [i] 
and y G [i + 1,?t.] because then, for some 01,02 G [m], M{y,ai) = M{y,a2) = 
1 — C{y) = 1 — by, M(l, oi) = 1 — M(l, 02) and we can see that K2 is isomorphic to 
M{{x, y}, {oi, 02, gy, hy}), which contradicts K2 ^ M({x, y}, ). So we have to assume 
that i <x <y <n. As M{[i + 1,/], 1) = C{[i + 1,1]), y > I. 



As K2 ^ M{{x, y}, ), no column of M{{x, y}, [m]) can equal C{{x, y}) = {b^., h 



\T. 



y> 



in particular, M{y, g^) = M{y, h^) = 1 — by (as M(x, g^) = M{x, h^) = b^). But then 

K2 = M{{l,y},{g,,h^,gy,hy}), (4) 
which is a contradiction proving our theorem. □ 

Theorem 13.21 yields that sat{n,K2) = forb(n, X2) =n+l which, in our opinion, 
is rather surprising. A greater surprise is yet to come as we are going to show now 
that sat(ri, K^) is constant for n > 4. 

Theorem 3.3 For K3 the following holds: 

7, zfn = 3, 



sat(n, K-:< 



10, ifn>A. 
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Proof. The claim is trivial for n = 3, so assume n > 4. A computer search [TO] 
revealed that 

sat(4, Ks) = sat(5, K3) = sat(6, Ks) = sat(7, K3) = 10, (5) 

which suggested that sat(?7,, X3) is constant. An example of a Ks-saturated 6 x 10- 
matrix is the following. 

"0000110111" 
1 1 1 1 1 

10 10 10 11 

1000011011 ■ 

1 1 1 1 1 
10 10 110 1 

It is possible (but very boring) to check by hand that M is indeed ii'3-saturated as is, in 
fact, any n x 10-matrix M' obtained from M by duplicating any row, cf. Theorem 12.11 
(The symmetries of M shorten the verification.) A /Cs-saturated 5 x 10-matrix can 
be obtained from M by deleting one row (any). For n = 4, we have to provide a 
special example: 

[00000001 1 1" 
1 1 1 1 1 

0011001101 ■ 

10 10 10 110 

So sat(n,iC3) < 10 for each n > 4 and, to prove the theorem, we have to show 
that no Ks-saturated matrix M with at most 9 columns and at least 4 rows can exist. 
Let us assume the contrary. 

Claim 1. Any row of M G SAT(n, ii'3) necessarily contains at least four O's and 
at least four I's, for n > 4. 

Suppose, on the contrary to the claim, that the first row M(l, ) contains only 3 
zeros, say in the first three columns. (By Lemma [3. II we must have at least 3 zeros.) 

If we replace the ith of these zeros by 1, i G [3], then the obtained column Cj, 
if added to M , does not create any K'3-submatrix. Indeed, the first row of [M, Cj] 
contains at most three O's while Cj([2,?T,]) is a column of M([2,?t,],). As M is K^- 
saturated, Ci, C2 and C3 are columns of M. These columns differ only in the first 
entry from M(,l), M(,2) and M(,3) respectively. Therefore, for each A G [2,?t,]'^'^\ 
the matrix M(v4, ) can contain at most e(M) — 3 < 6 distinct columns. But then 
any column C which is not the column of M and has the leading entry 1 (C exists 
as n > 4), contradicts the ii'3-saturation of M as the first row of [M, C] contains at 
most 3 zeros. This contradiction proves Claim 1. 

Therefore, e(M) is either 8 or 9. As we are free to complement the rows, we 
may assume that each row of M contains exactly four I's. Call A G [n]'-^-* (and also 
M(y4, )) nearly complete if M{A, ) has 7 distinct columns. 
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Claim 2. Any nearly complete M{A, ) contains (0, 0, 0)^ as a column. 

Indeed, otherwise M{A, ) D Tg-^ which already contains four I's in each row; this 
implies that the (one or two) remaining columns must contain zeros only. This implies 
that M{A, ) ^ K^, which is a contradiction. 

Claim 3. Every nearly complete M{A, ) contains as a submatrix. 

Indeed, if (0,0,1)^ is the missing column of M{A,), then some 7 columns span 
a copy of K3 \ (0, 0, 1)"^. By counting I's in the rows we deduce that the remaining 
column(s) must have exactly one non-zero entry and one of them equals (0, 0, 1)^, 
which is a contradiction. 

Now fix any nearly complete A which exists by the Ka-saturation of M. Assume 
that A = [3] and that the first 7 columns of M([3], ) are distinct. We know that the 
3-column missing from M([3], [7]) has at least two I's. 

If (1,1,1)"^ is missing, then M([3], [7]) contains exactly three ones in each row, 
so the remaining column(s) of M must contain an extra 1 in each row. As (1, 1, 1)"^ 
is the missing column, we conclude that e(M) = 9 and the 8th and 9th columns 
of M([3],) are, up to a row permutation, (0,0,1)"^ and (1,1,0)"^. This implies that 
M([3],) contains the column (0,0,0)"^ only once. On the other hand, by Claims 2 
and 3, M must contain the columns ((0)")"^ and ((0)'^"^, 1)"^ (as they cannot create 
K^) whose first three entries are zeros, which is a contradiction. 

Similarly, if (1, 1,0)^ is missing, one can deduce that, up to a row permutation, 
M([3],) contains 7 distinct columns plus the columns (1,0,0)-'" and (0,1,0)"^ and, 
again, the multiplicity of (0, 0, 0)"^ is only one, which is a contradiction as above, 
completing the proof of the theorem. □ 

We do not have any non-trivial results concerning Kk, k > 4 except that a com- 
puter search [10] showed that sat(5, K4) = 22 and sat(6, K4) < 24. (We do not know 
if a i^4-saturated 6 x 24-matrix discovered by a partial search is minimal.) 

Problem 3.4 For which k > 4, sat(n, Kk) = 0(1) ? 

4 Forbidding Small Matrices 

Here we will try to compute sat(n, F) for forbidden matrices with at most 3 rows. 
4.1 Forbidding 1-Row Matrices 

For any given 1-row matrix F, we can determine sat(n, F) for all but finitely many 
values of n: 
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Theorem 4.1 Let F = ((O)'", (1)') = [mT^, IT^] with I > m. Then, forn>l-l, 



Proof. Assume that / > 2, as the case / = 1 is trivial. 

Suppose that m = 0. The case / = 2 is trivial, so assume / > 3. An example of 
M G SAT(ri, F) can be built ioi n> l — lhj taking T° plus plus Xln]\{i}, i ^ [^~2], 
and Xli-2]- Clearly, each row of M D T° has exactly / — 1 ones so M is F-saturated 
and the upper bound follows. 

On the other hand, suppose on the contrary that some F-saturated matrix M has 
n > I — 1 rows and c < I columns. As c < 2" and M contains the all-0 column, c = I 
and some row M{i, ) contains at / — 1 ones. As we are not allowed multiple columns 
in M, some other row, say M{j, ), has at most / — 2 ones. Then X{j} is not a column 
of M because its ith entry is zero but its addition does not create / ones in a row. 
This contradiction establishes the case m = 0. 

For m > 1, let M consist of plus X{i}: ^ ["^ ~ 2], plus X[n]\{i}, i ^ " 1] 
and X[m-i,i-i]- Clearly, each row of M contains / ones and m — 1 zeros so any new 
column (which must contain at least one 0) creates an F-submatrix and the upper 
bound follows. The lower bound is trivial. □ 

Remark 4.2 The case when n < / — 2 in Theorem 14.11 seems messy so we do not 
investigate it here. 

4.2 Forbidding 2-Row Matrices 

Now let us consider some particular 2-row matrices. 

Let F = IT2, that is, F consists of the column (1, 1)-^ taken / times. Trivially, 
for / = 1 or 2, sat(n, /Tg^) = n + I with T^^ and [T^-'^,r"] being the only extremal 
matrices. For / > 3, we can only show the following lower bound which is almost 
sharp for / = 3 when we can build a 3r2^-saturated n x (2n + 2)-matrix by taking 
T^^, X[n-i], X[n], plus X{i,n} for i G [u - 1]. 

Lemma 4.3 For I > 3 and n>3, sat(n, IT^) >2n + l. 

Proof. Let M = [T^^,M'] be /ii'l-saturated. Note that M' must have the property 
that every column xa with A E [n]^'^^ either belongs to M' or creates an F-submatrix; 
in the latter case, exactly / — 1 columns of M' have ones in both positions of A. 
Therefore, adding to M' some columns of (possibly multiple), we can obtain a 
new matrix M" such that, for every A e [n]^'^\ M"(A,) contains the column (1, 1)^ 
exactly /—I times. If we let set Xj be encoded by the ith row of M" as its characteristic 




J + m 



1 , if m > 1 and I > 2. 



if I = 1 and n > 1, 
if m = and I = 2, 
if m = and I > 3, 
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vector, we have that |Xj fl Xj| = I — 1 for every 1 < i < j < n. The result of Bose [8] 
(see [T6l Theorem 14.6]), which can be viewed as an extension of the famous Fisher's 
inequality [12], asserts that the rows of M" are linearly independent over the reals 
or M" has two equal rows, say Xi = Xj. The last case is impossible because then 
= / — 1 and each other Xh contains Xi as a subset; this in turn implies that the 
column ((1)")-'" appears at least I — 1 > 2 times in M" and (since n > 3) the same 
number of times in M', a contradiction. Thus the rank of M" over the reals is n. 
Since every column added to M' when we were constructing M" was already present 
in M', the matrices M' and M" has the same rank over the reals. Thus M' has at 
least n columns and the lemma follows. □ 

Let us show that Lemma H73] is sharp for some n. Suppose there exists a symmetric 
{n, k,2)-design (meaning we have n fc-sets Xi, . . . ,Xn G [n]^''^ such that every pair 
{i,j} € [n]^'^^ is covered by exactly two Xj's). Let M be the n x n-matrix whose 
rows are the characteristic vectors of the sets Xi. Then [T^^,M] is a STg^-saturated 
n X (2n + l)-matrix. For n = 4, we can take all 3-subsets of [n]. For rz, = 7, we 
can take the family {[7] \ : z G [7]}, where Yi, . . . , I7 G [7]*-^'' form the Fano plane. 
Constructions of such designs for n = 16, 37, 56, and 79 can be found in [9], Table 
6.47]. 

Of course, the non-existence of a symmetric {n, k, 2) design does not imply any- 
thing about sat(n, 3T|), since a minimum 3T2^-saturated matrix [T^^,M] need not 
have the same number of ones in the rows of M. 

Lemma H73l is not always optimal for / = 3. One trivial example is n = 3. Another 
one is n = 5: 

Lemma 4.4 sat(5,3T|) = 12. 

Proof. Suppose on the contrary that we have a 3T2^-saturated 5 x 11-matrix [Tj-^, M]. 
Let Xi, . . . ,X^ be the subsets of [5] encoded by the rows of M. Then \Xi n Xj\ = 2 
for 1 < i < j < 5. If Xj C Xj for distinct i,j G [5], then \Xi\ = 2, every other 
Xh contains Xi as a subset, and M has two equal columns, a contradiction. In 
particular, 3 < \Xi\ < 4 for every i G [5]. A simple case analysis gives a contradiction 
by assuming that each \Xi\ = 3. Finally, if some \Xi\ = 4, say Xi = [4], then each 
of X2, . . . , X^ contains 5 and some two elements of [4], and we can easily derive a 
contradiction. □ 

Problem 4.5 Determine sat(n,3T2^) for every n. 

Remark 4.6 It is interesting to note that if we let F = [IT2, (0, 1)^] then the sat- 
function is bounded by a constant. Indeed, complete M' = [X[n]\{i}]i&[i] to an arbitrary 
F-saturated matrix M. Clearly, in any added column all entries after the Ith position 
are either zeros or ones; hence sat(n, F) < 2 ■ 2K 

It is easy to compute sat(T2,,T2^) by observing that the n-row matrix My whose 
columns encode Y C 2^^^ is T2^-free if and only if y is a chain, that is, of each two 



11 



members of Y one is a subset of the other. Thus My is Tg^-saturated if and only if Y 
is a maximal chain without repeated entries. As all maximal chains in 2["1 have size 
n + 1, we conclude that 

sat(n, T^) = forb(n, T^) =n+l, n>2. 
Theorem 4.7 Let F = [r°,r|] = 

Proof. For n > 3, the matrix M consisting of the columns (0, 1, (l)""^)"^, (1, 0, (1)""^)-^ 
and (0, 0, (l)""^)^ can be easily verified to be F-saturated and the upper bound fol- 
lows. 

Since n = 2 is trivial, let n > 3. Any 2-column F-free matrix M ^ T2 is, without 
loss of generality, the following: we have rows (0, 0), (1,1), (1, 0) and (0, 1) occurring, 
in this order, noo > 0, nu > 1, tt-iq < 1 and rioi < 1 times respectively. But then the 
addition of a new column ((0)"""^^, 1,1,...)"^ cannot create an F-submatrix. □ 

Theorem 4.8 Let F = = 

sat(n, F) = forb(n, F) = n + 1, n > 2. 

Proof. Clearly, forb(n, F) < forb(n, K2) = n + 1. 

Next, suppose that some M e SAT(n, F) has two equal rows, for example, 
M(l, ) = M(2, ) = ((1)', (O)'"). Let X = [I] and F = [/ + 1, / + m]. Define 

A = {j e[l + m]:M{i,j) = l}, ie[n]. 

(For example, Ai = A2 = X.) As M is F-free, for every i,j G [n], the sets Ai and 
Aj are either disjoint or one is a subset of the other. For i G [3,n], let fej = 1 if Ai 
strictly contains X or Y and let bi = otherwise (then Ai is contained in X ot Y). 
Let bi = 1 and 62 = 0. 

Clearly, C = (61, . . . , bn)'^ is not a column of M so its addition creates a forbidden 
submatrix, say F C [M, C]{{i,j}, ). Of course, bi = bj = is impossible as (0, 0)^ ^ 
F. If bi = bj = 1 then necessarily Ai n Aj ^ ^ and M{{iJ},) D (1, 1)^ contains F. 
Finally, if bi 7^ bj, e.g., bi = 0, bj = 1 and i < j, then Ai ^ Aj (as (0, 1)"^ cannot be a 
column of M{{i,j}, )), which implies Ai = Aj] but then we do not have a copy of F 
as (1, 0)"^ is missing. 

Thus, no F-saturated matrix M cannot have two equal rows and, by Theorem 12. 11 
sat(n, F) > 71 + 1 for any n. □ 

Trivially, sat(n, [(0, 1)"^, T2]) = 2 so we know the sat-function for any simple 2-row 
matrix. 



1 
1 



Then sat(ri, F) = 3, n > 2. 



1 1 

1 1 



Then 



12 



4.3 Forbidding 3- Row Matrices 



Here we consider some particular 3-row matrices. First we solve completely the case 
when F = [T0,T|]. 



Theorem 4.9 Let F = [T^,Ti] 



sat(n, F) 




Then 



3 or n > Q, 

4 or 5. 



Proof. For the upper bound we define the following family of matrices. 



Ma 



10 10 10 110 

110 110 10 

1 1 1 1 1 

1 1 1 

110 110 10 10 
10 110 10 110 
1110 110 1 
1 1 1 1 
0000000011 

10 10 10 
10 110 
110 10 
10 10 10 
1 1 1 
1 1 1 



[ ] for every 
e(M„) for 3 < n < 10. It 

(6) 



For any n > 7 let M„([6],) = Mg and M„(z,) 
7 < i < n. A computer search [10] showed that sat(?7,, F) 
remains to show that 

sat(n,F) > 7 

for n > 11. In order to see this, we show the following result first. 

Claim. If M is an F-saturated matrix of size n x m with n > 11 and m < 6 then 
M contains a row with all zero entries or with all one entries. 

Suppose on the contrary that we have a counterexample M. We may assume 
that the first 6 entries of the first column of M are equal to 0. Consider a matrix 
A = M([6], {2, . . . , m}). Note that every column of A contains at most two entries 
equal to 1, otherwise M([6], ) 3 F. Hence, the number of I's in A is at most 2(m — 1). 
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By our assumption, each row of A has at least one entry 1. Since 2(m — 1) < 12, 
A has a row with precisely one entry equal to 1. We may assume that A{1,1) = 1 
and A{l,i) = for 2 < i < m. Let C2 be the second column of M (remember that 



Consider the ri-column C3 = [0, C2({2, . . . , n})'^]'^ which is obtained from C2 by 
changing the first entry to 0. If it is not in M, then F C [M, C3]. This copy of F 
has to contain the entry in which C3 differs from C2. But the only non-zero entry in 
Row 1 is M(l, 2); thus F C [C2, C3], which is an obvious contradiction. Thus we may 
assume that C3 is the third column of M. 

We have to consider two cases. First, suppose that C2({2, . . . , n}) has at least one 
entry equal to 1. Without loss of generality, assume that 6*2(2) = (73(2) = 1. 

It follows that C2{i) = 6*3(2) = for 3 < i < 6 (otherwise the first and the second 
columns of M would contain F). Let 



By our assumption, each row of B has at least one 1; in particular m > 4. Clearly, 
B contains at most 2(m — 3) < 8 ones. Thus, by permuting Rows 3, ... ,6 and 
Columns 4, . . . , m, we can assume that .6(1, 1) = 1 while B{1, i) = for 2 < i < m— 3. 
Let C4 be the fourth column of M and C5 be such that 6*4 and C5 differ at the third 
position only, i.e., 6*4(3) = 1 and 6*5(3) = 0. As before, 65 must be in M, say it is 
the fifth column. Since 64({4, 5, 6}) has at most one 1, assume that 6*4(6) = 6*4(6) = 
6*5(5) = 65(6) = 0. We need another column Cq with Cq{5) = Cq{6) = 1 (otherwise 
the fifth or the sixth row of M would consist of all zero entries). In particular, m = 6. 
But now the new column 67 which differs from Cq at the fifth position only {i.e. 
6*7(5) = and Cj{i) = CQ^i) for i ^ 5) should be also in M, since M is F-saturated. 
This contradicts e(M) = 6. Thus the first case does not hold. 

In the second case, we have 62(2) = Cs{i) = for every 2 < i < 6. We may define 
i? as in ([7]) and get a contradiction in the same way as above. This proves the claim. 

Suppose on the contrary to the theorem that we can find an F-saturated matrix 
M with fi > 11 rows and m < 6 columns. By the claim, M has a constant row, say 
M'^ = [N"^ , (0)^™-*]"^. If 6* is an (n — l)-column missing from A^, then the column Q = 
{C^ , 1)^ is missing in M. Moreover, a copy of F in [M, Q] cannot use the n-th row. 
Thus F C [N, 6], which means that G SAT(ri - 1, F) and sat{n - 1, F) < m < 6. 
Repeating this argument, we eventually conclude that sat (10, F) < 6, a contradiction 
to the results of our computer search. The theorem is proved. □ 



6*2(1) =A(1,1) = 1). 



B = M({3,4,5 



6},{4,...,m}). 



(7) 



Theorem 4.10 Let F = [T^, r|, T|] 



111 
10 11 
110 1 



Then 
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Moreover, for any n > 8, sat(n, F) < 7. 



Proof. For n = 3, ... ,7 the statement follows from a computer search [10] with the 
following F-saturated matrices. 



Ma 



1 
1 



1 1 









1 1 





1 



1 1 
1 1 



1 1 1 1 



1 1 

1 










1 



1 





1 



1 

1 1 



111 
11 



1 1 1 1 



1 







1 1 



1 










1 



1 











11 

1111 



For any n > 7 let M„([6],) = Mq and M„(z,) 
7 < i < n {i.e. the last row of Mg is repeated (n 
M, 

contains no pair of equal rows. 
Conjecture 4.11 Let F = [r3°,r|,T|]. 



:[ 1 1 1 1 ] for every 
6) times). It remains to show that 



„, n > 8, is F-saturated. Clearly, this is the case, since Mj is F-saturated and F 



□ 



Then sat(n, F) = 7 for every n > 8. 



Theorem 4.12 Let F = 



<2 












1 





Then 



sat(n, F) 



7, ifn = 3, 
10, z/4<n<6. 



Moreover, for any n > 7, sat(n, F) < 10. 

Proof. For n = 3, . . . , 6 the statement follows from a computer search [TU] with the 
following F-saturated matrices. 



Ma 



10 10 110 1 

1 1 1 1 1 

1 1 1 1 1 

1 1 1 



15 



1 1 1 1 1 

10 110 10 1 

1 1 1 1 1 

1 1 1 1 1 

1 1 1 1 1 



For any n > 6 let M„([5], ) = Mg and M„(i, )=[ 1 1 1 1 1 ] for 
every 6 < i < n. It remains to show that M„, n > 7, is F-saturated. Clearly, this is 
the case, since Mg is F-saturated and F contains no pair of equal rows. 



□ 



Conjecture 4.13 Let F = Tg- . Then sat(n, F) = 10 for every n > 7. 



Theorem 4.14 Let Fi 



Then sat(n, Fi) 
sat(n, Fi) < 3n 




1 
1 



andF2 = [T|,r|] 




1 
1 



= sat(?7,, F2) = 3n — 2 for any 3 < n < 6. Moreover, for any n > 7, 
2 and sat(n, F2) < 3n — 2 as well. 



Proof. Let M„ = [r°, T^,T^, T^], where C and consists of all those columns of 
which have precisely one entry equals 1 either in the first or in the nth row, e.g., 
for n = 5 we obtain 



Ms 



0100001111000 
0010001100100 
0001001010010 
0000101001001 
0000011000111 



Clearly, e(M„) = e{T^) + e{T^) + e{T^) + e{f^) = l+n + l + 2n-A = 3n-2. Moreover, 
since is Fi-admissible we get that M„ is both Fi and F2 admissible. Now we show 
that Mn is Fi-saturated (consequently, M„ is also F2-saturated) . Indeed, pick any 
column C which is not present in M„. Such a column must contain at least 2 ones 
and 1 zero. Let l<i,j,k<nhe the indices of C so that q = 0, Cj = = 1. If z = 1 
or i = n, then the matrix M„({z, j, fc}, ) contains Fi. Otherwise, Ci = c„ = 1, and 
there also exists 1 < i < n such that Cj = 0. Here M„({1, i,n},) contains Fi. Hence, 
sat(n, Fi) < 3n — 2 and sat(n, F2) < 3n — 2 for any n > 3. A computer search [TU] 
yields that these inequalities are equalities when n = 3, . . . , 6. □ 



Conjecture 4.15 Let Fi 

3n — 2 for every n > 7. 



Tl and F2 = [Fg^Ff]. Then sat(n,Fi) = sat(n,F2) 



Remark 4.16 It is not hard to see that sat{n, Fi) > n + c^fn for some absolute 
constant c and all n > 3. Indeed, let M be an n x (n + 2 + A) Fi-saturated matrix of 
size sat(n, Fi) for some A = \{n). We may assume that M(, [n + 2]) = [T°, F^, T"]. 



Moreover, we assume that every column of matrix M([A],{n + 3, 



n 



2 + A}) 
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contains at least one entry equal to 1 (there must be a permutation of the rows 
of M satisfying this requirement). We claim that all rows of M({A + 1, . . . , n}, {n + 
3, . . . , n + 2 + A}) are different. Suppose not. Then, there are indices A + 1 < z, j < 
such that M(i, {n + 3, . . . , n + 2 + A}) = M(j, {n + 3, . . . , n + 2 + A}). Now consider a 
column C in which the only nonzero entries correspond to i and j. Clearly, C is not 
present in M, since the first A entries of C equal 0. Moreover, since M is Fi-saturated, 
the matrix [M, C] contains Fi. In other words, there are three rows in M which form 
Fi as a submatrix. Note that the and jth row must be among them. But this is 
not possible since Fi has no pair of equal rows. 

Let Mo = M({A + 1, . . . , n}, {n + 3, . . . ,n + 2 + A})^. Clearly, Mq is Fi admissible. 
Anstee and Sali showed (see Theorem 1.3 in [3j) that forb(A, Fi) = O(A^). That 
means that n — \ = O(A^), and consequently, A = Vl{^/n). Hence, sat(n, Fi) = 
e(M) >n + Vt{^/n), as required. 

References 

[1] R. P. Anstee, Forbidden configurations: Induction and linear algebra, Europ. J. 
Combin. 16 (1995), 427-438. 

[2] R. P. Anstee, J. R. Griggs, and A. Sali, Small forbidden configurations, Graphs 
Combin. 13 (1997), 97-118. 

[3] R. P. Anstee and A. Sali, Small forbidden configurations IV: The 3 rowed case, 
Combinatorica 25 (2005), no. 5, 503-518. 

[4] R.P. Anstee and Z. Fiiredi, Forbidden submatrices. Discrete Math. 62 (1986), 
225-243. 

[5] B. BoUobas, Combinatorics, Set Systems, Families of Vectors, and Combinatorial 
Probability, Cambridge Univ. Press, 1986. 

[6] B. BoUobas, Extremal graph theory. Handbook of Combinatorics (R. L. Graham, 
M. Grotschel, and L. Lovasz, eds.), Elsevier Science B.V., 1995, pp. 1231-1292. 

[7] J. A. Bondy, Induced subsets, J. Combin. Theory (B) 12 (1972), 201-202. 

[8] R. C. Bose, A note on Fisher's inequality for balanced incomplete block designs, 
Ann. Math. Statist. 20 (1949), 619-620. 

[9] C. J. Colbourn and J. H. Dinitz (eds.). Handbook of combinatorial designs, sec- 
ond ed.. Discrete Mathematics and its Applications (Boca Raton), Chapman & 
Hall/CRC, Boca Raton, PL, 2007. 

[10] A. Dudek, O. Pikhurko, and A. Thomason, On minimum saturated matrices (it 
includes the C source code of satmat . c), e-print arXiv:math/xxxxxxxx. 



17 



[11] p. Erdos, A. Hajnal, and J. W. Moon, A problem in graph theory, Amer. Math. 
Monthly 71 (1964), 1107-1110. 

[12] R. A. Fisher, An examination of the possible different solutions of a problem in 
incomplete blocks, Ann. Eugenics (London) 10 (1940), 52-75. 

[13] P. Frankl, On the trace of finite sets, J. Combin. Theory (A) 34 (1983), 41-45. 

[14] P. Frankl, Z. Fiiredi, and J. Pach, Bounding one-way differences. Graphs Combin. 
3 (1987), 341-347. 

[15] Z. Fiiredi, Turdn type problems. Surveys in Combinatorics, London Math. Soc. 
Lecture Notes Ser., vol. 166, Cambridge Univ. Press, 1991, pp. 253-300. 

[16] S. Jukna, Extremal combinatorics with applications to computer science. Springer 
Verlag, 2001. 

[17] O. Pikhurko, The minimum size of saturated hypergraphs, Combin. Prob. Com- 
puting 8 (1999), 483-492. 

[18] N. Sauer, On the density of families of sets, J. Combin. Theory (A) 13 (1973), 
145-147. 

[19] S. Shelah, A combinatorial problem: Stability and order for models and theories 
in infinitary languages, Pac. J. Math 4 (1972), 247-261. 

[20] A. Sidorenko, What we know and what we do not know about Turdn numbers. 
Graphs Combin. 11 (1995), 179-199. 

[21] V. N. Vapnik and A. Chervonenkis, The uniform convergence of frequences of the 
appearance of events to their probabilities (in Russian), Teor. Veroyatn. Primen. 
16 (1971), 264-279. 



18 



C source code of satmat 



* Copyright (C) 2009 Andrzej Dudek, Oleg Pikhurko and Andrew Thomason 
* 

* This progrcun is free software: you can redistribute it and/or modify 

* it under the terms of the GNU General Public License as published by 

* the Free Software Foundation, either version 3 of the License, or 

* (at your option) any later version. 
* 

* This progrcun is distributed in the hope that it will be useful, 

* but WITHOUT ANY WARRANTY; without even the implied warranty of 

* MERCHANTABILITY or FITNESS FDR A PARTICULAR PURPOSE. See the 

* GNU General Public License for more details. 
* 

* You should have received a copy of the GNU General Public License 

* along with this prograim. If not, see <http://www.gnu.org/licenses/>. 
* 

* 

* satmat. c: this progrsun finds the smallest number of columns (saturated 

* matrix) with n rows without having an r-row matrix F 
* 

* compile: 

* °L gcc -Wall -ansi -04 -o satmat satmat. c 
* 

* compile w/ debug option: 

* 7, gcc -Wall -ansi -04 -o satmat -DDEBUG satmat. c 



* 



usage : 



* 



7, . / satmat 



* 



* 



* Authors : Andrzej Dudek 

* Oleg Pikhurko 

* Andrew Thomason 



adudekOandrew . emu . edu 

pikhurkoOandrew . emu . edu 

A . G . ThomasonQdpmms .CEim.ac.uk 



* 



*/ 



#include <stdlib.h> 
#include <stdio.h> 
#include <time.h> 
#include <string.h> 



#define MAXN 10 



/* max value of n */ 



#define MAXR 5 



/* max value of r */ 



#define MAXNCR 252 /* should be MAXN choose MAXR */ 



#define MAXS 50 



/* max size of strong sat set */ 
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int sub_mask[l << MAXN] [MAXNCR] ; /* the sub-column mask indexed by some r rows, 

e.g. n=5 r=3 subset 4 in colex is 10011 
30th column is 11110 subcolumn = 110 = 6 
so subcol [30] [4] = ..111011111 */ 

int subs_left [MAXS] [MAXNCR] ; /* allowed sub-columns indexed by rows found so 

far; e.g. if the first 18 columns under 
consideration contain, in rows indexed by 
subset 4, sub-columns of types 0,1,4,5 
then subs_left [18] [4] = "110011 = "51 */ 

int poss_cols[(l « MAXN) * MAXS]; /* list of columns to try next */ 

int T[MAXR] [1 « MAXR] ; /* family of T_k'l matrices */ 

int F[MAXR] [1 << MAXR] ; /* forbidden matrix */ 

int colT [MAXR+2] ; /* colT[l] is the first column of T_k~l in T */ 

int irrel [1 « MAXN] ; /* no sub-column is in F */ 

int mask; /* column representation of matrix F */ 



void gen_T(int r) 
{ 

int rows [MAXR] ; 
register int i, 1, col; 

/* init T */ 

memset((void *)T, 0, MAXR*(1 « MAXR) *sizeof (int) ) ; 

colT[0] = 0; 
col = 1; 

for (i = 0; i < r; i++) -[ /* precisely (i+1) Is per column */ 

colT[i+l] = col; 

for (1=0; 1 < i+1; 1++) { 
rows [1] = 1; 

} 

do { 

for (1=0; 1 < i+1; 1++) { 
T[rows[l]] [col] = 1; 

> 

col++; 

/* find next subset of i Is */ 
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for (1=0; 1 < i; 1++) { 

if (rows[l] + 1 != rows[l + 1]) 
break; 

> 

> 

rows [1] += 1 ; 

while (—1 >= 0) { 
rows [1] = 1; 

> 

} while (rows[i] < r) ; 

} 

colT[i+l] = col; 

> 

void gen_mask(int r, int c) 
i 

register int i, j, col; 

#ifdef DEBUG 

printf("\n F mask\n"); 
#endif 

mask = 0; 

for (j = 0; j < c; j++) i 

for (col =0, i = 0; i < r; i++) { 
if (F[i] [j] == 1) { 
col 1= (l«i) ; 

} 

} 

mask I = (1 << col) ; 
#ifdef DEBUG 

printfC mask[7.d] — > 7.d\n", j, col); 

#eiidif 
} 

#ifdef DEBUG 

printfC mask — > 7od\n" , mask); 
#endif 
} 

void gen_irrel (int n, int minO, int mini, int *i 
{ 

register int i,j,ones; 

#ifdef DEBUG 

printfC \n minO — > 7.d\n" , minO) ; 

printfC mini — > 7.d\n" , mini); 
#endif 
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/* column with less than mini Is or more than (n-minO) Is */ 
for (*irr =0, i = 0; i < (1 « n) ; i++) { 

for (ones = 0, j = 0; j < n; j++) { /* count Is */ 
if (i & (1 « j)) { 
ones++ ; 

> 

} 

if (ones < mini I I ones > n-minO) { 
irrel [(*irr)++] = i; 

} 

} 

#ifdef DEBUG 

printf("\n irrelevant columns — > "); 
for (i = 0; i < *irr; i++) { 
printf("7.2d ", irrel [i] ) ; 

} 

printf ("\n") ; 
#endif 
} 

int is_irrel(int col, int irr) 
{ 

register int i; 

/* binary search could be more effective */ 
for (i = 0; i < irr; i++) i 
if (irrel [i] == col) { 
return 1 ; 

} 

> 

return 0; 

} 

void input (int *n, int *target, int *r, int *c, int *irr) 
{ 

int i , j ,k, col ,minO , mini ; 
char ans ; 

printf (" \n************************** Wellcome to SATMAT! **************************\n") ; 
printf ("This program checks if sat(n,F) <= target, for a given an r-row matrix F.\n"); 
printf ("If so, an n-row F-saturated matrix is produced. \n") ; 

printf ("\n input n — > "); 

scanf ("7.d" , n) ; 

if (*n < 2 II *n > MAIN) { 

printf (" n is not in the range 2 to %d (MAXN)\n\n", MAXN) ; 

exit(O) ; 
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> 



printf("\n input r — > "); 

scanf ("7.d" , r) ; 

if (*r < 2 II *r > MAXR) { 

printfC r is not in the range 2 to 7A (MAXR)\n\n", MAXR); 

exit(O) ; 

} 

if (*r > *n) { 

printfC r is not in the range 2 to 7A \n\n", *n) ; 
exit(O) ; 



printf("\n target — > "); 
scanf ("7.d" , target) ; 

if (*target < 2 II *target >= MAXS) { 

printfC target is not in the range 2 to "/.d (MAXS-1) \n\n" , MAXS - 1) ; 
exit(O) ; 

> 

printf("\n input F as a union of Ti matrices\n\n") ; 
gen_T(*r) ; 

for (i = 0; i < *r; i++) { 

for (j = 0; j <= *r; j++) { 
if (i == (*r-l)/2) { 

printfC T7,d=|", j); 

> 

else { 

printfC I"); 

> 

for (k = colT[j]; k < colT[j+l]; k++) { 
printf Cy.d", T[i] [k]) ; 

} 

printfC I ") ; 

} 

printf ("\n") ; 

> 

col = 0; 

minO = mini = *r; 

for (i = 0; i <= *r; i++) { 

printf ("\n add T7,d (y/n) — > ", i) ; 

scanf (" 7oC", &ans) ; 

if (ans == 'y') { /* copy Ti into F */ 
for (j = 0; j < *r; j++) { 

for (k = colT[i]; k < colT[i+l]; k++) { 
F[j] [col+k-colT[i]] = T[j][k]; 

} 

> 

col += (colT[i+l]-colT[i] ) ; 
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/* at least this number of Os in very column*/ 
if (*r-i < minO) { 
minO = *r-i; 

> 

/* at least this number of Is in very column*/ 
if (i < mini) { 
mini = i; 

> 

> 

} 

*c = col; 

/* print F */ 
printf ("\n") ; 

for (i = 0; i < *r; i++) { 
if (i == (*r-l)/2) { 
printf (" F=|"); 

} 

else { 

printf (" I"); 

> 

for (j = 0; j < col; j++) { 
printf ("%d", F[i] [j]); 

} 

printf (" I \n") ; 

} 

/* find F mask */ 
gen_mask(*r, *c) ; 

/* find irrelevan columns */ 
gen_irrel (*n, minO, mini, irr) ; 

} 

void setupCint n, int r, int c, int *ncr) 
{ 

int rows [MAXR] ; 

register int 1, col, subr, subcol; 

/* subset of r rows */ 
for (1=0; 1 < r; 1++) { 
rows [1] = 1; 

> 

/* index of subset in colex */ 
subr = ; 
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do { 

/* for each column, extract subcolumn indexed by rows*/ 
for (col = 0; col < (1 « n) ; col++) { 

for (subcol =0, 1=0; l<r; 1++) { 
if (col & (1 « rows[l])) { 
subcol 1= (1 « 1) ; 

} 

> 

sub_mask[col] [subr] = (1 « (1 « r)) - 1 - (1 « subcol); 
sub_mask [col] [subr] &= mask; 

> 

/* find next colex subset of rows */ 
for (l=0;l<r-l; 1++) { 

if (rows[l] + 1 != rows[l + 1] ) { 
break; 

> 

> 

rows [1] += 1 ; 
while (—1 >= 0) { 
rows [1] = 1; 

} 

subr += 1 ; 

} while (rows[r - 1] < n) ; /* don't do colex beyond n */ 

*ncr = subr; /* we stopped at subset n choose r */ 
f f lush(stdout) ; 

} 

void print_sol (int *cols, int n, int 1, int irr) 
{ 

int matrix [MAXN] [MAXS] ; 
int i , j ; 

/* columns from the solution */ 
for (j = 0; j < 1; j++) ■[ 

for (i = n-1; i >= 0; i— ) { 

matrix [i] [j] = ((l«i) & cols[j]) ? 1 : 0; 

} 

> 

/* irrelevant columns */ 
for (j = 0; j < irr; j++) { 

for (i = n-1; i >= 0; i— ) { 

matrix [i] [1+j] = ((l«i) & irrel[j]) ? 1 : 0; 

> 

> 

1 += j; 
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for (i = 0; i < n; i++) { 



if (i == (n-l)/2) { 



printfC M_{y.d,7.2d>=| ", n, 1); 

} 

else { 



printfC I"); 

} 



for (j = 0; j < 1; j++) { 

printfC "7.d" , matrix[i] [j]); 

} 

printfC I \ii") ; 

} 

f f lush(stdout) ; 

> 

void find_sat (register int n, register int r, register int ncr, register int target, int irr) 



register int 1, m, subr, col; 
register int *p, *q, *s, *t , *u; 

register int strong = target +1; /* smallest strong set so far */ 
register int examples = 0; 

int gone_critical [MAXNCR] ; /* sub-columns which became critical */ 
int *gc [MAXS] ; /* pointers into gone_critical */ 

int crit_at [MAXNCR] ; /* when a sub-column went critical */ 

time_t start, stop; 
double diff; 
#ifdef DEBUG 

long long counter = 0; 



{ 



int cols [MAXS] ; 

int *ptr [MAXS] , *to [MAXS] ; 



/* current set of columns */ 
/* pointers into poss_cols */ 



#endif 



start 



time (NULL) ; 



for (m 
if 



= 0, col = 0; col < (1 « n) ; col++) { 
(is_irrel (col , irr) == 0) { 
poss_cols [m++] = col; 



} 



> 



ptr[0] 
to[0] 



= poss_cols; 
• poss_cols + m; 



#ifdef DEBUG 

printf("\n number of relevant columns — > %d\n" , m) ; 
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f f lush(stdout) ; 
#endif 



/* there is only one relevant column */ 
if (m == 1) { 

print_sol((int*)0, n, 0, irr) ; 

strong = 0; 

examples = 1; 

> 

gc [0] = gone_critical ; 
for (m = 0; m < ncr; m++) { 
crit_at [m] = -1; 

> 

1 = 0; 

/* we have set of 1 + 1 columns */ 
do { /* main loop */ 

#ifdef DEBUG 

counter++; 

#endif 

cols[l] = *ptr[l]++; 

if (1 + 1 >= strong I I ptr[l] > to[l]) i 

1 -= 1; 

if (1 < 0) { 
break; 

> 

goto undo; 

} 

#ifdef DEBUG 

/* info */ 
if (1 < 5) { 

for (m = 0; m <= 1; m++) { 
printf ("7.3x" , cols[m]); 

> 

put char ( ' \r ' ) ; 
f f lush(stdout) ; 

} 

#endif 

/* make up subs_left mask after column 1 */ 
if (1 == 0) { 

for (subr = 0; subr < ncr; subr++) { 

subs_lef t [0] [subr] = sub_mask [cols [1] ] [subr] ; 

> 

p = subs_lef t [0] ; 
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} 

else { 

p = subs_left[l - 1]; 

} 

q = subs_left[l] ; 

s = sub_mask [cols [1] ] ; 

t = p + ncr ; 

u = gc[l] ; 

while (p < t) { 

m = *q++ = *p++ & *s++; 

if ((m & (m - 1)) ==0) { /* only bad subcolumn */ 
m = ncr - 1 - (t - p) ; 

if (crit_at[m] < 0) { /* subcolumn not critical so far */ 
crit_at [m] = 1; 
*u++ = m; 

} 

> 

} 

gc[l + 1] = u; 

p = ptr[l]; /* find columns which extend this set */ 
q = s = to [1] ; 
t = gc[l + 1] ; 
while (p < s) { 
m = *p++; 

u = gc[l]; /* check only subcolumns which went critical */ 
while (u < t) { 
subr = *u++; 

if ((subs_left [1] [subr] & sub_mask [m] [subr] ) == 0) { 
m = 0; 
breeik; 

} 

> 

if (m) { 

*q++ = m; 

> 

} 

ptr [1 + 1] = to[l] ; 
to[l + 1] = q; 

/* if possible, go round again with larger 1 */ 
if (ptr[l +1] < to[l + 1]) { 

1 += 1; 

continue ; 

} 

/* if not, set is inextendible upwards; is it strongly saturated? 
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subr = ; 



for (col=0; col < cols [0] ; col++) { 
if (is_irrel (col , irr) == 1) { 
continue ; 

> 

for (subr = 0; subr < ncr; subr++) { 

if ((subs_left [1] [subr] & sub_mask [col] [subr] ) == 0) { 
break; 

} 

> 

if (subr == ncr) { 
break; 

} 

} 

if (subr < ncr) { 

for (m = 0; m < 1; m++) { 

for (col = cols[m] + 1; col < cols[m + 1]; col++) { 
if (is_irrel (col , irr) == 1) -[ 
continue ; 

> 

for (subr = 0; subr < ncr; subr++) { 

if ((subs_left[l] [subr] & sub_mask [col] [subr] ) == 0) { 
break; 

} 

} 

if (subr == ncr) { 
break; 

} 

} 

if (subr == ncr) { 
break; 

> 

> 

} 

if (col == cols[l]) { /* we have a strong set size 1 + 1*/ 
if (strong < 1 + 1) { 

continue; /* not best possible */ 

> 

if (strong > 1 + 1) { 

strong =1+1; /* best so far */ 
examples = 0; 

> 

examples += 1; 

if (exEunples == 1) { 
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priiit_sol(cols , n, strong, irr) ; 

stop = time (NULL); 

diff = diff time (stop, start); 

printfC'Xn after — > 7.. If min\n", dif f /60) ; 

ff lush (stdout) ; 

> 

} 

/* now go again with ssune 1 */ 
undo : 

u = gc [1] ; /* sub-cols critical at 1 no longer are */ 
t = gc[l + 1] ; 
while (u < t) { 

crit_at [*u++] = -1; 

> 

} while (1 >= 0) ; /* main loop */ 
#ifdef DEBUG 

printf("\n number of iterations — > 7olld\n" , counter); 
#endif 

stop = time (NULL); 

diff = difftime(stop, start); 

printf("\n time taken — > 7,. If min\n" , dif f /60) ; 

printf ( " \n******************************* Solution ********************************\n") ; 
if (examples > 0) { 

printf ("7.28s sat(%d,F) = 7.d\n" , "",n, strong+irr) ; 

} 

else { 

printf ("7.28s sat(7.d,F) > 7.d\n" , "", n, target+irr) ; 

} 

} 

int main(int argc, char *argv[]) 
{ 

int n, r, c, ncr, target, irr; 

input (&n, fetarget, &r, &c, &irr) ; 

setup (n, r, c, fencr) ; 

find_sat(n, r, ncr, target-irr, irr); 

exit(O) ; 
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